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Abstract 

We prove a global existence theorem for the 3x3 system of relativistic Euler 
equations in one spacial dimension. It is shown that in the ultra-relativistic limit, 
there is a family of equations of state that satisfy the second law of thermodynamics 
for which solutions exist globally. With this limit and equation of state, which 
includes equations of state for both an ideal gas and one dominated by radiation, 
the relativistic Euler equations can be analyzed by a Nishida-type method leading to 
a large data existence theorem, including the entropy and particle number evolution, 
using a Glimm scheme. Our analysis uses the fact that the equations of state are 
of the form p = p{n,S), but whose form simplifies to p = a^p when viewed as a 
function of p alone. 
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CHAPTER 1 

Introduction 

1. The Compressible Euler Equations 

The compressible Euler equations form a nonlinear system of first order partial 
differential equations that models continuous medium. Nearly seventy 

years after Newton wrote down the laws of motion in his Principia for a system 
of discrete particles, F = ma, Euler and d'Alembert produced a linear, continuum 
theory of sound waves. These sound waves obeyed the linear wave equation, 

□m = Utt — (?Dw{u) = 0, 

where c > is the sound speed. Several years later, Euler wrote down the evolution 
equations for the nonlinear theory of sound waves. Today these equations are written 
as 

Pt + Div [pu] = 0, 
{pu)t + Div [pu ® M + pi] = 0, 
(1) Et + Dtv [{E + p)u] = 0, 

where subscripts in the independent variables denotes partial differentiation and 
Div = d/dx + d/dy + d/dz. In three spacial dimensions, the compressible Euler 
equations ([1]), also called Euler's equations, form a system of five equations with 
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§ L 1 . THE COMPRESSIBLE EULER EQ UATIONS 

Figure 1.1. The "breaking" of a wave front which produces a shock wave. 

six unknowns, p, e, u\ and p, which closes when an equation of state, p = p{p, S), 
is prescribed. In the following we will focus our study on the case of one spacial 
dimension. Under this assumption, the Euler equations reduce to a system of three 
equations: 

Pt + Hx = 0' 

ipu)t+[pu^+p]^ = 0, 
(2) E, + [{E + p)ul = 0. 

It is well known that even for smooth initial data, discontinuities form in the fluid 
variables in the solution to the Cauchy problem in finite time, [3]. Qualitatively, the 
nonlinearities in the equations cause waves to propagate at different speeds leading 
to the "breaking" of waves. See Figure 11.11 This loss of regularity corresponds to 
the emergence of shock waves. 

The Euler equations are a particular example of a system of conservation laws. 
A system of conservation laws in one spacial dimension is a first order quasi-linear 
system of partial differential equations of the form 




(3) 



Ut + F{U), = 0, 



§1.1. THE COMPRESSIBLE EULER EQ UATIONS 

where U = (f/i, . . . , f/„) are the conserved quantities and F{U) = (Fi(f/), . . . , Fn{U)) 
the fluxes. Much of the early work on the general structure of systems of conservation 
laws was set out by Lax, \5\. Lax's results provided the foundation necessary for 
Glimm to give the flrst general existence theorem in 1965, |4]. Glimm's fundamental 
result provided a new way to analyze shock wave interactions. In the 1990's, Bressan, 
Liu and Yang headed a push for the well posedness of the general n x n Cauchy 
problem, [2]. 

A Nishida system is a speciflc class of conservation laws, which in certain cases 
includes the the Euler and Relativistic Euler equations, that allows one to prove 
global existence of solutions. In particular, the shock-rarefaction curves in a Nishida 
system behave nicely in the large. Nishida and Smoller were flrst to gave a global, 
large initial data, existence proof for the compressible Euler equations with a particu- 
lar equation of state, [8]. Shortly after this. Temple extended Nishida and SmoUer's 
global existence result by including the entropy evolution of the gas, [12]. More 
recently, Smoller and Temple proved that under certain conditions the Relativistic 
Euler equations also form a Nishida system, |10] . 

It should be noted that the existence theorem for a general system of conservation 
laws comes at a cost; we require the initial data to be of sufficiently small total 
variation. The smallness requirement is needed because the structure of the shock- 
rarefaction curves can exhibit complicated nonlinear phenomenon in the large. When 
sufficiently small data is considered, the analysis can be confined within a small region 
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§1.2. THE RELATIVISTIC EULER EQUATIONS 

in state space in which the shock-rarefaction curves have a canonical structure that 
can be exploited when analyzing solutions. 



2. The Relativistic Euler Equations 

In 1905, Einstein introduced the special theory of relativity. Within this frame- 
work one can generalize the classical Euler equations to obtain equations that fit 
within the theory of relativity. 

The relativistic compressible Euler equations in one spatial dimension form a 
system of three equations, 

{u''n),a = 0, 
(4) = 0, (3^0,1, 

where T"''^ is the stress energy tensor for a perfect fluid. 



and the subscript ", a" denotes partial differentiation with respect to the coordinate 
x". We will use Einstein's summation convention where repeated up-down indices 



§1.2. THE RELATIVISTIC EULER EQ UATIONS 
are summed and adopt the following notation: 

Components of the 2- Velocity 

p Proper Rest Energy Density 

p Pressure 

e Specific Internal Energy 

n Bar yon Number 

S Specific Entropy 

T Temperature 

The components of the Minkowski metric 77"^ are given by 



-1 



For convenience, we will also use units where the speed of light is unit, c = 1. 
The proper energy density, p, is related to the particle number density and the 
internal energy by p = n{l + e), [14J. This equation is the sum of the rest mass 
energy nc^ = n and the internal energy ne. Furthermore, thermodynamics provides 
a functional relationship between the quantities, e, T, S, p and n. This relationship 
is given by the second law of thermodynamics, [3]: 



(5) de = TdS + ^dn. 
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§1.2. THE RELATIVISTIC EULER EQUATIONS 

The relativistic Euler equations (jlj) can be written as a system of conservation 
laws by choosing a particular Lorentz frame and writing the instantaneous worldline 
trajectory of the fluid, m", in terms of the classical velocity v. The components of 
(m^jM^) are proportional to the vector (l,f) and is of unit length according to the 
inner-product defined by the metric rj. From this we find the components are 
related to v by 

Using this, the first equation is equivalent to 

d f n \ d f nv \ 

The second and third equations in (jl]) can also be rewritten. With /5 = we find 
T^^a = gives 

and with (3 = 1, T^^^, = gives 

Simplifying the terms inside, we can write the system (jlj) as the system of conserva- 
tion laws. 



(6) 



Ut + F{U)^ = 0, 



§1.2. THE RELATIVISTIC EULER EQ UATIONS 
where, 



c = iTr^p +p' ^-^ 



and 



(8) F{U)=(-==,{p + p)- -+p,{p + p) 

It is interesting to note that the relativistic Euler equations are indeed a gener- 
aUzation of the classical Newtonian equations of hydrodynamics ([2]). To see this we 
view (j6]) under the assumptions of a classical fluid; fluid velocities are small com- 
pared to the speed of light and the pressure is dominated by the rest mass. More 
specifically, we assume |f | ^ 1 and p/p <^ 1. Under these assumptions the equations 
T°'^a = become the equations of motion of a classical gas: 



9 d , , ^ 



and 

|(P«) + |(P«^+P)=0. 

Notice that the density of the fluid in the classical Euler equations is now replaced 
by the proper mass-energy density. The new variable n is used for conservation of 
particle number. 

Nearly all terrestrial phenomenon falls into the classical, Newtonian case. In 
a hurricane, for example, wind speeds may reach speeds of dOm/s. However, this 
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§1.2. THE RELATIVISTIC EULER EQUATIONS 
velocity is insignificant when compared to the speed of light, 

\v\ ~ 90m/s ~ 3 X 10"^c = 10"^ < 1. 

Furthermore, the pressure to mass density ratio, p/ p, can be shown to be of the order 
of 10^^^, \Jj. In this situation the classical Euler equations would certainly suffice. 

It is clear from the last example that even in seemingly extreme situations on 
Earth, they are far from relativistic events. We must look to the cosmos to find 
examples where a gas has a high enough pressure to make p/p non- negligible and 
sufficiently high velocity to make the relativistic correction terms such as — 
important to the gas' evolution. These situations arise in astrophysical events such 
as gamma-ray bursts, solar flares and in remnants of supernovas. The relativistic 
Euler equations are also used in modeling the early universe, |13] . 

Like the classical Euler equations, the relativistic Euler equations are not closed; 
an equation of state relating thermodynamic variables is needed to close the system. 
This choice of equation of state changes the characteristics of the evolution of the 
gas and has a significant effect on the complexity of its analysis. A natural equation 
of state for a gas is one satisfying the ideal gas law and whose internal energy is 
proportional to its temperature. Using the second law of thermodynamics, one finds 
the relation 

(9) e(n,5) = e^V-\ 

which for some constant 7 > 1 is called a polytropic equation of state. A polytropic 
equation of state is typically used to model air in the classical sense with 7 ^ 1.4. It 
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1.2. THE RELATIVISTIC EULER EQUATIONS 



is known that using this equation of state vacuums may form in a solution to ([2]) and 
(jl]) when velocities and densities are sufficiently large to completely void a region of 
matter. Vacuums pose problems in the standard estimating techniques and at this 
point prevents one from obtaining large data existence theorems, [9]. 

A class of equations of state one typically encounters which still include most 
desirable dynamics are called barotropic, given hj p = p{p). The class of equations 
of state, p = a? , for 1 < 7 < 2 are barotropic and are used in astrophysical 
modeling. In this case < a is constant, [T]. 

If one considers (11 01) . the limiting case of the barotropic equation of state p = a'^p"' 
when 7 = 1, the system ([6]) contains special properties; in this limit one can prove 
global solutions exist for initial data with arbitrarily large, but finite, total variation, 
[10] . Moreover, vacuums do not form in the solution. 



In this thesis we will extend these results to prove large data existence theorem 
for an ultra-relativistic gas with an equation of states of the form 



(10) 




(11) 



e{n,S) = A{S)n 



7-1 



where the function A satisfies the following conditions: 



(Al) 



A : M" 



+ 



(A2) 



A e C^(M+), 
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§1.2. THE RELATIVISTIC EULER EQUATIONS 
(A3) A'{S) > for 5 > 0. 

The family flTTl) includes equations of state for a polytropic gas and one dominated 
by radiation satisfying 



12 e{n,S) = 

n 



Using the relation p = n{l + e), the equations of state (fTT!) do not reduce to (fTOj) . 
However, they do in the ultra- relativistic limit. For the ultra-relativistic limit, we 
assume the internal energy dominates the rest mass energy; in other words, p = ne. 
Under this assumption, an equation of state of the form ffTTj) reduces to an equation 
of state of the form ffTOl) with = (7 — 1). We take advantage of the fact that 
in this limit the pressure is still a function of n and S, but whose form reduces to 
(JTOl) when viewed as a function of p alone. This model now allows one to find the 
entropy and particle number density evolution of the gas and still take advantage of 
the simphfying effects of an equation of state of the form (fTOl) . 

The particular equation of state (fT2l) is also used to model massless thermal 
radiation. In this case the ultra-relativistic assumption is not needed since the mass- 
energy in p drops out, leaving only the internal energy. In particular for 7 = 4/3, the 
radiation dominated equation of state is used to model the early universe, because 
this radiation has been predicted to make the dominant energy contribution, |13j . 
In either situation, massless particles or in the ultra-relativistic limit, we still have 
an equation of state of the form (ITUl) . 
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It is interesting that for the classical Euler equations there is only one way to 
assign an entropy profile to a gas with an equation of state of the form (|TOl) . This 
equation of state is given by 



(13) e{p,S) = aHn{p) + ^ + C, 



for constants a > 0, C > 0. A global existence theorem for the classical Euler 
equations with this equation of state was given by Temple in 



3. Statement of Main Theorem 

The goal of this paper is to prove the following: 

1. Theorem. Let po{x), Vo{x) and So{x) be arbitrary initial data satisfying, 
Po{x) > 0, -1 < fo(x) < 1 andSoix) > 0. Let S = \n[A{S)] for e{n, S) = A{S)n"'-\ 
1 < 7 < 2, and A satisfying (lAip . (]A2p and (lA3p . Suppose further that 

(14) VariEoi-)} < oo, 

(15) \/ar{ln(po(-))} < oo, 
and 



(16) Var{ In ( ] } < oo. 



l-Vo{-) 
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§1.3. STATEMENT OF MAIN THEOREM 

Then there exists a bounded weak solution {p{x, t),v{x, t), S{x, t)} to in the Ultra- 
Relativistic limit, satisfying 

(17) \/ar{S(-,t)}<iV, 



(18) Var{\n{p{-,t))} < N, 

and 

where N is a constant depending only on the initial variation bounds in (1141) . f|T5l) . 
and (fT6|). 

It should be noted that Theorem [1] is a generahzation of the work by SmoUer and 
Temple in |10] that includes the entropy evolution. In other words, in this model 
we are able to prove global solutions exist including a physically relevant entropy 
and particle number density profile. SmoUer and Temple found that the relativistic 
Euler equations with equation of state (fTO!) possessed the property that after each 
elementary wave interaction in a Glimm scheme, Var{ln{p)} is non-increasing. This 
functional, introduced by Liu, is used as a replacement for the quadratic potential in 
Glimm's original analysis, which can be used to show that f|T5|) and f|T6|) implies f|T8l) 
and f|T9l) . Considering the ultra-relativistic limit, the solutions of Riemann problems 
are independent of the value of S, enabling one to solve for the intermediate state 
in the projected state space and place a corresponding entropy wave between them. 



13 

§L3. STATEMENT OF MAIN THEOREM 

In |10] it is shown that for an equation of state of the form (|TOl) . the shock curves 
are translationally invariant in the plane of Riemann invariants. In our case, this 
property continues to liold under certain coordinate changes in the three dimensional 
non-projected state space for an equation of state of the form (ITT]) . This can be 
viewed as the relativistic analogue of the large data existence result in |12] with a 
family of distinct entropy profiles. 

The main part of the analysis is showing that Var{S} is bounded in our ap- 
proximate solutions. We extend the analysis by SmoUer and Temple for the ultra- 
relativistic regime with equation of state given by f|TT]) . by utilizing the geometry 
of the shock curves in the space of Riemann invariants. If we only considered the 
variation of S across shock waves, we find that Var{S} is uniformly bounded by 
Var{ln{p)} for a polytropic equation of state. However, across the linearly degen- 
erate entropy waves, there is no change in pressure, and hence no jump in proper 
energy density by ffTOl) . Thus, another method must be employed to estimate the 
strengths of these jumps. For a gas dominated by radiation or for a general equation 
of state of the form fill I) , the situation seems more dire as the change in entropy 
across a shock depends on the initial entropy value. It is not known a priori that 
this dependence does not lead to blow-up in the variation in S. 

Furthermore, in certain elementary wave interactions, Var{S} may actually in- 
crease while Var{ln{p)} remains invariant. Complicating matters, using Aln(p) as 
the definition of wave strengths increases the technicality of the entropy wave es- 
timates. For example, after the interaction of two shocks of the same family, the 
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entropy change across the new shock may be less than the sum of the jumps across 
the two preceding shock waves. This happens because the new shock wave has 
strength less than the sum of the two previous. In other words when two shock 
waves combine, the strengths are not simply additive, but the new wave strength is 
strictly less than the simple sum of the incoming shock strengths. It follows that 
under certain circumstances the change in entropy across the new single shock may 
be less than the sum of the entropy jumps across the approaching shocks. 

To alleviate these technicalities, we propose a more classical approach by using 
the change of Riemann invariants as a measure of wave strength. More specifically, 
the strength of a 1(3)— shock is determined by the change in the first (third)-Riemann 
invariant and the contact discontinuity by the change in entropy or a specific change 
in a function of the entropy. Using the change in Riemann invariants as a measure 
of wave strength for a Nishida system was used to prove existence of solutions in [6] , 
[8] and [12] . Under this regime, wave strengths are now additive and the sum of all 
the strengths of shock waves is shown to be non-increasing in time. Moreover, the 
wave interaction estimates can be analyzed as in the classical case. In conclusion, 
using A ln(p) as a measure of wave strength dramatically simplifies the interaction 
estimates for the nonlinear waves, but complicates the problem dealing with entropy. 

In summary, we show there exists a family of equations of state, which include 
the case of a polytropic and radiation dominated gas, that one can use and obtain 
a global existence theorem. These equations of state allow one to also calculate 
the entropy and particle density associated with the gas. This is in contrast with 
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the classical case where there is one equation of state with the same properties 
corresponding to very heavy molecules. 

The rest of this paper is outlined as follows: 

In Chapter [21 we give a detailed analysis of the structure of simple wave solutions 
of Q). Using these properties, we prove global existence of solutions to Riemann 
problems. Furthermore, we obtain a priori wave interaction estimates which will 
be used to produce estimates on approximate solutions constructed using a Glimm 
scheme in Chapter [31 

In Chapter [31 we give an overview of the Glimm difference scheme and prove 
estimates on the approximate solutions obtained for system ([6]). Chapter [H contains 
the proof of our main theorem. 
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CHAPTER 2 
Relativistic Gas Dynamics 
1. Gas Dynamics 

We consider a gas wliere the proper energy density and pressure satisfy the 
relationship (ITO!) . Causahty restricts the sound speed Cg = \/ d-pjdp = a to be 
less than unity. Under assumption (fTOj) . the system ([6]) decouples so that we may 
solve for two variables first, then solve for the third afterward. In this section, we will 
show in the domain p>0, — l<f<l, and 5 > 0, Riemann problems are globally 
solvable and their general structure consists of two waves separated by a jump in 
entropy traveling with the fluid. We then discuss wave interaction estimates which 
will allow us to prove global existence of solutions using a Glimm scheme in Chapter 
m Our analysis uses the special geometry of the shock and rarefaction curves in the 
space of Riemann invariants. 

To begin, we will compute the eigenvalues and eigenvectors associated with the 
relativistic Euler equations ([6]). In order to simplify this process, we will exchange 
the first equation, conservation of particle number, with the equivalent equation that 
says that entropy is constant along flow lines. We note that this equation holds for 
continuous solutions, but fails when shock waves form since entropy increases across 
shocks, [3]. 
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2.1. Proposition. For smooth solutions of (Q, the following supplemental equa- 
tion holds: 

(20) m"S',„ = 0. 

More specifically after choosing a particular Lorentz frame, 

(21) St + vS, = Q. 

Proof. We show that conservation of energy and momentum is equivalent to 
continuous flow being adiabatic, i.e. fl2Tl) . We take the stress-energy tensor of a 
perfect fluid, 



T"^ = n (l + e + ^) M V + pr/"^, 
= nuju^u^ + pr]"^^ , 

with uj = (l + e + ^) for convenience. Then conservation of energy-momentum 
equation, T"^^ = 0, is given by 



n - 

= {nuu''u^)^^+p,f,r'', 
(22) = nn^(^n")^^+p,^r/"^. 
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§2.1. GAS DYNAMICS 

where conservation of particle number, (^nu^^ ^ = 0, is used in the last step. Multi- 
plying (l22l) by —Ua and summing we find, 

(23) = —UaT'^^p = ~nu^ {uu") f^Ua— P,l3U^- 

To simplify this expression, we claim. 

Indeed, 

Ua {OJU"") ^ = UaU°'0J^p + (MqU"^) UJ = —UJ^p, 

where the second term, (m^m"^^) u, vanishes because 

= = 2moM°^. 



Thus, ( l23l) now reads. 



= nuj^pu^ — p^pu^ 



n I e^pu^ + ( v^pu^ + V u^ ) - p^^jw^. 



,/3 



/3y 



nTS^pu^ . 



The last step follows from the second law of thermodynamics. Since n, T 7^ we 
conclude, u^S^p = 0. Furthermore, after choosing a particular frame of reference and 
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replacing the worldline trajectory with 



1 



u = 



we get 



1 




V 



S. 



y/1 — v' 



x 



0. 



In particular, since 



1 



7^ 0, dsn) holds. 



□ 



It is interesting to note that (l20l) continues to hold in curved spacetimes within 
general relativity. For this case, differentiation is replaced by covariant differentia- 
tion. 

In order to solve the Riemann problem by a series of simple waves, we need to 
know that the corresponding wave speeds are distinct. If this is the case the system 
is strictly hyperbolic. 

2.2. Definition. We call a system of conservation laws ([3]) Strictly Hyper- 
bolic in an open connected subset U MP if at each point u E U , dF has n real 
distinct eigenvalues, {Aj(u)}"^^, such that 



Since a strictly hyperbolic system has n distinct eigenvalues, the corresponding 
eigenvectors form a basis at every point in U . Along with strict hyperbolicity, we 
require one more assumption on the eigenvector-eigenvalue pairs; the corresponding 
eigenvalues are either constant or monotonically increasing or decreasing along the 
integral curves determined by the eigenvectors. 



Ai(n) < . . . < \n{u). 
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2.3. Definition. Let {{^i{u), Ri{u))}^^^ be the eigenvalue-eigenvector pairs as- 
sociated with dF for a strictly hyperbolic conservation law in an open connected subset 
?7 C M" with \i{u) < . . . < \n{u). We call the characteristic field Genuinely 
Non-Linear in U if for all u E U, 

R,{u) ■ VXiiu) ^ 0, 

and Linearly Degenerate if for all u eU , 

Ri{u) ■ V\i{u) = 0. 

In the following proposition we characterize the three eigenclasses of the system 

2.4. Proposition. Let p = a^p with < a < 1. Then the system ^ is strictly 
hyperbolic at (p, f , S) for p > —1 < v < 1 and 5 > 0. Furthermore, the first and 
third characteristic fields are genuinely non-linear and the second linearly degenerate. 

Proof. Equivalent systems of equations possess the same eigenvalues, so we will 
replace the conservation of particle number equation with the equivalent equation 
(12T]) . Since the flux functions ([8]) are complicated implicit functions of the conserved 
variables ([7]), our plan is to rewrite the conservation laws ([6]) as 



LOt + G{uj)uj^ = 0, 
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where lo = (p, v, S)'^, then calculate the eigenvalues and eigenvectors in terms of 
these variables. To do this we rewrite ([6]) using the chain rule as 

A{Lj)ujt + B{ijj)ujx = 0, 

then find G{uj) by multiplying on the left by A~^ to get 



+ [A-'^B] {uj)uj, = 0. 



By the chain rule, 



A{co) 



2v 



(1-1)2)2 



and 



V 

(«' + 1)t^ ia' + l)pj^ 
Note that A{uj) is invertible because for — 1 < t> < 1 and p > 0, 

Det[A(a;)] = ^ V 7^ 0. 



(1 



After some work we get 



A-\uj) 









2v 



1+1)2 



1— a2D 



(l-i;2)(l+a2i;2 „3_ 



(a2+l)p{l-a2t;2) p(l-a2D2) 
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Therefore, 



G{uj) = [A-^B] 



.00 



(a^+l)P Q 
a2(i-t,2)2 {a?-l)v „ 



We look for the roots of the characteristic polynomial, 



= Det [G{uj) - A/] 



(v - A)(A(-1 + av)-a + v){-a - v + \{1 + av)) 
{av — 1)(1 + av) 



There are three values of A that make the numerator zero. 



(24) 



Ai 



V — a 
1 — fa' 



Ao = V, As 



V + a 
1 + va 



We show for < a < 1 and — 1 < f < 1, 



Ai < A2 < A3 



Indeed, 



v"^ < 1 <^==^ —av'^ > —a <^==^ V — av^ > v — a. 



By the restrictions on v and a, (1 — av) > and thus, 



f (1 — av) > V — a <^=^ V > 



V — a 
1 — av 



Showing V < {v + a)/{l + va) is similar, we omit the details. We conclude that for 
p>0, — l<f<l, and S > the system is strictly hyperbolic. 

Now, we show that the first and third characteristic fields are genuinely nonlinear 
and the second is linearly degenerate. To do this we need to find the eigenvectors of 
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G{uj). For A2 we simply find 

R2{p,v,S) = (0,0,1)' 

and after some work, 



and 

^a(l-w2)^ 

Computing tlie gradients of tlie eigen- fields with respect to u; = (p, v, S), 

VA2 = (0,1,0), 
Thus, 

(1 — Of )^ 
i?2-VA2 = 0, 

i?3-VA3 = ~ ^ 0. 

(1 + av)^ 

The first and third are non-zero and bounded by the restrictions on a and v. □ 



It is interesting to note that the eigenvalues fl2^ are the relativistic analog of the 
sum of the local sound speed and fluid velocity in the classical Euler equations. In 
the classical case, the first and third characteristic fields have eigenvalues Ai = u — c 
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and A3 = u + c, which are the sum and differences of the ffuid speed and the local 
speed of sound respectively. The eigenvalues fl2^ are exactly the relativistic sum of 
two velocities within the frame work of relativity. 

2. Riemann Invariants 

The Riemann invariants for the system ([6]) can be found from the eigenvectors. 
An i**^— Riemann invariant is a function ijj such that 



In other words, the level curves of ip are the integral curves of the i characteristic 
field. We will perform our interaction estimate analysis in the coordinate system of 
Riemann invariants because the rarefaction curves have particularly simple structure; 
straight lines parallel to the coordinate axes. From the eigenvector Ri we see that 
along 1— rarefaction curves, 



i?. . VV' = 0. 



dp 
dv 



+ 1 p 



a 1 — v' 



9 ' 
1^ 



which we can explicitly solve to find that along the first integral curve. 
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This can be done for the third integral curve in a similar fashion. We therefore 
define: 



The function r = r(p, v) is constant across 3— rarefaction waves and s = s{p, v) is 
constant across 1— rarefaction waves. From the supplemental equation fl2Tl) . we see 
that the entropy, S*, is a third Riemann invariant constant across 1 and 3— rarefaction 
waves. In our analysis, we will view state space in the coordinates of the Riemann 
invariants rather than the conserved variables. However, using S is not sufficient 
because the shock curves in (r, s, S) space are, in general, not translationally invari- 
ant. Instead we will use E = \\i{A{S)) as our third coordinate. It will be shown in 
Section H] that in (r, s, S) space, the shock-rarefaction curves are indeed independent 
of base point. Since S" is a Riemann invariant, \n{A{S)) must be one too. Indeed, 
suppose that ip is a. i**^— Riemann invariant and let / G C"'^(]R, M). Then f{ip) is an 
i*^— Riemann invariant as well since. 



We now change our variables from the conserved quantities (f/i, ?72, f/3) to (p, f, 5*). 

2.5. Proposition. In the region, p > 0, —1 < v < 1, S > 0, the mapping 
{p,v,S) — > {Ui, 1/2,113) is one-to-one, and the Jacobian determinant of the map is 
both continuous and non-zero. 




(25) 




Ri ■ V/(^) = fWRi ■ VV^ = 0. 
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Proof. We will show first that the map (p, v) {U2, U^) is one-to-one for p > 
and — 1 < i; < 1. Assume the contrary. Suppose we have (pi,fi) and {p2,V2) such 
that U2{pi,vi) = U2{p2,V2) and U^^pijVi) = Us{p2,V2)- To begin we show that if 
Vl = V2 = V then pi = p2. From the equality U^i^pi^v) = Uz{p2,v) we have 

Pi ((a^ + 1)^ + l)=P2 ((a^ + l)y^ + 1) . 

Since the term 

((a^ + 1)^ + 1)^0 

for any — 1 < 1; < 1 we must have pi = p2- We now show that if the images of U2 and 
U3 are equal, then we must have vi = V2 and, by the previous argument, pi = p2- 
From U3{pi,vi) = U^{p2,V2) and U2{pi,vi) = U2{p2,V2) we have 

(26, 

and 

Pi^f Vl \ f V2 \ 
P2 \l-vlJ \l-viJ- 

Note that if fi/(l — vf) = we must also have f 1 = and by fl26l) . f 2 = since 

P1/P2 7^ 0. Assume that Vi 7^ 0. 

Replacing pi/p2 in ( l26l) and simplifying, 

a'^vf + 1 _ a^t;! + 1 

Vl V2 

which further reduces to 

(fi - V2){a'^ViV2 - 1) = 0. 
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Since |a|, \v2\ < 1, the second term, (a^t>it>2 — 1) 7^ 0, so it must be vi = V2- 
Therefore, the mapping (p, f) < — > (f/2, U3) is one-to-one. 

Now we show that the mapping {p,v,S) {Ui,U2,Uz) is one-to-one. Proceed 
again by contradiction by supposing (pi, f 1, Si) and (p2, "^2, 5'2) have the same image. 
Since U2 and f/3 only depend on p and f , the previous argument shows that pi = 
P2 and Vi = V2. We now show that Si = 5*2. Since n = n{p, S) the equahty 
f/i(pi,fi, 5*1) = f/i(p2,f2, 5*2) reduces to 

n(p, Si) = n(p,5'2). 

Therefore, we are done if dn/dS 7^ 0. We use the fact that p = ne to rewrite the 
second law of thermodynamics ([5]) as 

ndp = n^TdS + (a^ + l)pdn. 



Therefore, 



dn 
dS 



7^0, 



(a2 + l)p 

and the mapping (p, v, S) — > ([/i, U2Us) is one-to-one. 
The jacobian matrix of the map is given by 
/ 



J 



-I 



(1-^2)3/2 
V (l+a2)p 



a2 + i)£il±!^ 



a2 + l)^ + l 







whose determinant is 



det(J) 



n^T(lj-aV) 
(l-t;2)2 



>0, 
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which is continuous onp>0, — l<f<l and S > 0. 

□ 



3. Jump Conditions 

Systems of conservation laws, or more specifically the relativistic Euler equations 
(ini), encode the required information to calculate the evolution of discontinuities, i.e. 
shock waves, in one or more of the conserved variables. One must use care however, 
because systems of equations equivalent to for smooth solutions can, and typically 
do not, give the same relations for discontinuous solutions. A prime example of this 
is specific entropy is constant along fiow lines for continuous solutions of ([6]) from 
(I2T]) . but entropy is not conserved and increases across a shock front. 

For systems of conservation laws, the relations defining the dynamics of shock 
waves are the Rankine-Hugoniot jump conditions. These relations state for a shock 
wave traveling at speed s, the change in the conserved quantities U across the shock 
and the change in F{U) across the shock, denoted [[U]] and [[F(f/)]] respectively, 
satisfy, 

(27) s[[U]] = mU)]]. 

For a given state Ul, the Rankine-Hugoniot relations, for each i = 1, . . . ,n, define 
a 1— parameter family of states that can be connected on the right by a shock wave 
in the i*^ characteristic family. Moreover, this curve has second order contact with 
the curve defining all the states that connect to Ul on the right by an rarefaction 
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wave given by the i^^ integral curve. These facts were first proven by Lax in 1957 for 
a general system of strictly hyperbolic conservation laws with genuinely nonlinear or 
linearly degenerate characteristic fields, [5]. 

We call IZiiU) the integral curve of the i^^ characteristic field that passes through 
the state U and Si(U) the one parameter family of states defined by (127|) that defines 
states that connect to t/ by a shock wave in the i*^ family. Only half of each of these 
curves will be physically relevant. For a genuinely non-linear characteristic field we 
take the portion of TZiiJJ) extending from U that satisfies XiiJJ) < \i{U'). Call this 
portion IZf(U). On the other hand, take the portion of the shock curve SiiU) that 
satisfies the Lax entropy condition, 

UU') < s < \{U). 

Call this portion S^{U). Finally, define TiiU) = njiU) U S;{U). 

For our system given by ([6]), we have that the tangent to the worldline of the 
shock front is proportional to (1, s). Define by 

(/°,/^) = (!,.). 

The jump conditions ( |271) for the system (jl]) is then given by 

[[nw"]] la = 0, 

(28) [[T°T]L = 0, P = 0,1. 
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Recall that 1^ is found by contracting with the metric rj: 

From the first equation in ( l28l) we have for some constant m, 
(29) m = nu^la = niulla- 

For the case m = 0, we have for n, > 0, 



U la — UlJa- 



Since the components are in a one-to-one relation with the fluid velocity we 
have V = v^. Furthermore, the second equation reduces to p = p^. This case, m = 0, 
corresponds to an entropy wave rather than a compressive shock. Shock waves will 
correspond to m 7^ 0. The thermodynamic relationships across a shock wave in a 
solution to the relativistic Euler equations was first given by Taub, 



2.6. Proposition (Taub, 1948). Let U = {p,v,n) and Ul = {pl.vl.ul) he two 
states separated by a shock wave. Then the following relation holds: 



/qn^ P + Pf < \ PL + PL I , ^ 
30 — Y- P + Pl) = 2 — (pL + P) ■ 



Proof. We will show that across a shock wave the following condition holds on 
the two separating states, 

(31) (p±iy.(pj^)\(p,-pUp^+E±±iA=o, 

\ n J \ riL J \ n ul J 
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Assuming this holds, we multiply out, cancel and collect terms with n and ni in 
the denominator on the left and right respectively to get 

p'^ + PP + ppl + PPl _ pI + PlPl + PPl + PlP 

Equation ( l30l) follows directly. 
For convenience define 

P + P -, Pl + PL 

9 = and = • 

n riL 

If m = we have a jump discontinuity. Since we are concerned about the shock 
waves, assume m 7^ 0. In this case the second equation in (!28|) gives 

ngu°'u^la + pr?°^/a = nLgLU^u^la + PlV'^^L, 
that, in light of (l29l) . reduces to 

(32) mgu'^ + pf = mgLu1 + pLl'^. 

Contracting equation (132|) with and then using (!29l) we find 

, P , PL 

(33) - g + - = gLUlUfS + — 

and 



(34) gu^UL, + ^ = -gL + —- 

riL riL 



n n 



We use ( l33l) to solve for m^M/?: 
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Plugging (|35|) into fl34|) . combining and using the definition of (7 and ql, we obtain 



Tlie global structure of the solutions of the shock relations fl7r|) for the relativistic 
Euler equations in the space of Riemann invariants was first done by SmoUer and 
Temple for an equation of state of the form (fTOj) . [10] . We summarize their results 
in the following lemma: 

2.7. Lemma (Smoller, Temple, 1993). Letp = a'^p withO < a < 1. The projection 
of the i-shock curves for i = 1,3 onto the plane of Riemann invariants (r, s) at any 
entropy level satisfy the following: 

(1) The shock speed s is monotonically increasing or decreasing along the shock 
curve Si and for each state {pl,vl) 7^ {priVr) on Sj the Lax entropy condi- 
tion holds: 



(2) The shock curves, when parameterized by Aln(p), are translationally invari- 
ant. Furthermore the 1 and 3— shock curves based at a common point (r,s) 
have mirror symmetry across the line r = s through the point (f , s) . 

(3) The i— shock curves are convex and 



□ 



In particular, with p = a'^p, (l30ll reduces to 



(36) 




MPR,Vr) < Si < Xi{pL,VL)- 
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for i = 1 and 




< 1 



ds -^2K-l 
for i = 3 where K = 2a^/ (1 + a^)^. 

In light of Lemma [?7n we see that we can globally define the shock curves Si{U) in 
the rs— plane and we know that everywhere on this curve the Lax entropy conditions 
hold. We now extend the analysis of SmoUer and Temple and show that the entropy 
change along the shock waves also possess the translationally invariant property 
and are convex in a particular coordinate system. After this we will show that 
the Riemann problem is globally solvable with equation of state (ITTl) . in the ultra- 
relativistic limit. 



In this section, we will show certain properties hold for our family of equations 
of state. Namely, we will need that as a function of wave strength, the change in a 
certain function of entropy is independent of base point. Moreover, we will find that 
the change of this function of entropy and its derivative are monotone increasing. 
We will use these facts in our estimates on the entropy waves in Section [61 

For an equation of state of the form 



4. Equations of State 



e(n, S) = A{S)n^-\ 
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with A satisfying flAl|) . flA2|) and flASp . the second law of thermodynamics says, 



2<9e 



p(n, S) = n — = (7 — l)y4(S')n''' = (7 — l)en. 
on 



In the ultra-relativistic hmit this further reduces to 



p{n, ^) = (7 - l)p, 



an equation of state of the form ( fTOl) with a = 1/7 — 1. 

Now we will show that a certain function of entropy across a shock wave is 
independent of base state {pl,vl, Sl) by using Proposition 12. 6[ Choose S by 

(37) S(5) = ln(A(5)). 

Our goal is to show that across a shock wave, the difference [S — S^] is a function 
of the change of the corresponding Riemann invariants alone. Then the difference 
[S — Sl] along the shock curve is independent of base point. Finally, we will show 
that the difference [S — S^] and its derivative, as a function of the change of Riemann 
invariants, are monotone increasing. Later we will measure the strength of 1— shocks 
as the change in r and by the change in s for 3— shocks. It is sufficient to show 
that the change [S — S^] and its derivative are monotone increasing as viewed as a 
function of ln(p/p2.), because they satisfy the relationship as parameters, 

2a ^ , , . 
Ar = ^— -Alnp . 

+ 1 



For 3— Shocks we replace Ar with As. Thus, 



d[S - Sl] _ d[S - Sl] 
d{r-rL) dln{p/pL) 



d\n{p/pi 



d{r - tl) 



c? + 1 d\S- Sl] 



2a d\n{p/pL) 
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Using fl36l) . we calculate [S — S^,], 

(38) S-Si = \n{AiS))-\niA{SL)), 



n-y PL 



5 



:i-7)ln(^) - Jin 



PL J 2'"U + (7-l)^ 



Thus for cr = ln(p/pi), 

(39) [S - Si](a) = (1 - 7)a + J In ^ ^ + ~ ^^^'^ 



2 Vl + (7-l)e"'' 
After differentiating, we have 

...^ - Sl] ^ (e-- 1)^(2 -7)(7-l) 

^ ^ da 2(l + e-(7-l))(e-+(7-l))' 

which is non-negative in the domain 1 < 7 < 2 and a > 0. Furthermore, the 
derivative is zero only when a = 0. Thus, [S — SL]{a) is a monotone increasing 
function. 

Next we show that d[T, — T,L]/da is also monotone increasing. We take another 
derivative and find 

d^[S-S^] _ ^2(2-7)(7-l)(e3--e-) 
da^ 2(1 + e-(7 - l))2(e- + (7 - 1))2 ' 

The denominator is always positive and the numerator is positive because 1 < 7 < 2 
and e^'^ > e'^ for a > 0. 
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We have proven the following proposition: 

2.8. Proposition. Consider the ultra-relativistic Euler equations with the equa- 
tion of state e{n, S) = A{S)n'^^^ and A satisfying (lAip . (]A2|) and (]A3|) . Then the 
change in T, = \n{A{S)), when regarded as a function of the change in the corre- 
sponding Riemann invariant, is independent of base state. Geometrically, the shock 
curves, as viewed in {r, s,!!)— space, are translationally invariant. 

An interesting fact is that the change in S becomes nearly linear for strong shock 
waves. We state this corollary. 

2.9. Corollary. Under the assumptions of Proposition \2.8[ the change in S 
becomes nearly linear for strong shocks. 

Proof. This follows immediately when considering the following limit: 

^.^ ^[S - S,] ^ (e-- 1)^(2 -7)(7-l) ^ (2-7) 

fT^oo da fT^oo 2(1 + e'^(7 - l))(e'^ + (7 - 1)) 2 

□ 

In the following sections we will show that both an ideal gas and one dominated 
by radiation fall into the family of equations of state given by (ITT!) . 



4.1. Ideal Gas. An ideal gas satisfies the ideal gas law, 

P=RT, 
n 
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[3]. Furthermore, if we assume that the internal energy is proportional to the tem- 
perature, 

-T, 



7-r 

we can use the second law of thermodynamics to determine e(n, S). From 

— = T = ^ ~ ^ 
dS 

we get the for some function ip, 

e(n, S) = e^^"^e 



To find if{n) we again use the second law of thermodynamics to get the relation 

de 7 — 1 



dn n 



-e, 



which reduces to 



dip 7 — 1 
dn n 



Solving for ip, we get 



ip{n) = (7 — 1) In (n) = In {n' 



Therefore, 



e(n, S) = e^V-\ 



for some 1 < 7 < 2. Thus, in the case of a polytropic gas, we have an equation of 
state of the form e{n, S) = A{S)n'^~^, where 



A{S) =e^^. 
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satisfying 1^ and (jAS]). We also see that S is proportional to S: 



= In(eV^) = ^S. 

Notice that in the case of an ultra-relativistic polytropic gas, it would have been suf- 
ficient to consider (r, s, S*)— space since the shock curves would still be translationally 
invariant. 

4.2. Radiation Dominated Gas. A gas in local thermodynamical equilibrium 
with radiation when only the internal energy and pressure are dominated by radiation 
is characterized by 



(41) e = and p = ^anT^, 

n 6 



where = 7.56 x 10 is the Stefan-Boltzmann constant, pLj. We can generalize 
the equation of state ( 14T1) to the continuum of equations of state. 



(42) e=— and p = h - l)aRT~ , 

n 

for 1 < 7 < 2. Notice (H2|) reduces to fj4T|) when 7 is chosen to be 4/3. 

In order to find the entropy profile associated with this equation of state we again 
use the second law of thermodynamics. From de/dn = p/n^ we find 

which in light of P2|) and after some algebra, reduces to 
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Similarly, from de/dS = T, we find 

{^)annTi^)§ 



71^ 



which can be simplified and integrated to find that for some function /(n), 
(44) -fanT^^) =nS + f{n). 

Differentiating (jH]) with respect to n and using ( l43i) . we find the following relation 
on /, 

nf'in) = fin). 

Thus, for some constant c, f{n) = cn, and we can incorporate c into the entropy 
level S giving, 

7a/{T^~) = nS, 

or equivalently, 

n 

Therefore, in the ultra-relativistic regime or for a massless gas. 



S 



P = Or I I rC 



and 



(45) e{n,S) = aR(^\ n^-\ 

The equation of state for thermal radiation ( l45l) is of the form ( ITT]) with 

A{S) =an[ 

\lCiR 
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Unlike the case for a polytropic gas, where it would have been sufficient to con- 
sider just the change in S, the change in entropy across a shock wave is no longer a 
function of ln{p/pL) alone; it is also dependent on the starting entropy level. More 
specifically, in this case we have the ratio of entropy values being independent of 
base state, rather than the difference, leading to 



Choosing the new coordinate S = ln(y4(S')) is necessary to keep the change indepen- 
dent of base point. In the case of a radiation dominated gas, 



Riemann problems are used as the building blocks of finite volume method solu- 
tion schemes for systems of conservation laws. The Riemann problem is a particular 
class of Cauchy problems with initial data of the form. 



We will show that for any to initial states in the region p > 0, — 1 < f < 1 and S > 0, 
there exists a solution of the Riemann problem for the system ([6]) with equation of 
state (ITTl). 





5. The Riemann Problem 




X > 0. 
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2. Theorem. Consider left and right states Ul = {pl,vl, Sl) andUR = {pr,vr, Sji), 
such that Pl,Pr > 0, — 1 < vi,vr < 1, and Sl,Sr > 0. With the equation of state 
( fTTl) satisfying 1 < 7 < 2, ( lAll) . (lA2p and ( ]A3I) . ^/iere exists a weak solution to 
the Riemann problem < Ul,Ur > for system in the ultra-relativistic limit. This 
solution is unique in the class of solutions with constant states separated by centered 
rarefaction, shock and contact waves. 

Proof. For any entropy level, the projection of the shock-rarefaction curves onto 
the rs— plane is translationally invariant by Lemma [2.71 We will show first that for 
any two states, Ul = {pl,vl) and Ur= {pr,vr) in the rs— plane, there exists a 
intermediate state Um = {pa4,vm) such that Um is on the shock-rarefaction curve 
based at Ul and Ur is on the shock-rarefaction curve based at Um- For convenience, 
let Ti{U) denote the projection of the i*^— shock-rarefaction curve based at U at any 
value of S onto the rs— plane. Given a state U l, partition the rs— plane into four 
regions: /, consisting of all states above Ti{U l) and to the right of %{Ul)', II, states 
above Ti{Ul) and to the left of T^iJJ l)', HI, states below Ti{U l) and above T^{U l)', 
and IV , states below Ti{Ul) and to the left of %{Ul)- See Figure [2TT1 

Consider Ti{U l)- For each U m G Ti{Ul) the 3-wave rarefaction curves based at 
Um extend vertically upwards, parallel to the s— axis. Therefore, for all states Ur 
in region / or //, there is a unique state Um G T^iiU l) that connects f/^ to fZ/j by 
a 1— shock or a 1— rarefaction wave followed by a 3— rarefaction wave. 

We now turn our attention to the portion below Ti(U l) in the rs— plane. For 
region IV, we notice for any state f/i G TVI{Ul) the shock curve S^{Ui) is a 
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II 



I 




Figure 2.1. A partition of the rs— plane into four sections: /, //, 
/// and IV. 

horizontal translation of S^{Ul)- Thus, all the 3— shock curves extending from 
1Z'^{U l) cover region IV . For region III it is clear that the 3— shock curves extending 
downward from {Ul) must cover all states in the region. But, we must show that 
if we take two states Ui and U2 on the shock curve of U l they will never intersect. 
Suppose two shock curves intersect at a third state f/3. See Figure [2^ We know by 
Lemma 12.71 that 



z V2ir - 1 

- < , — < 1. 

y -V2K - 1 
However, if Ui and U2 are on the same shock curve, 

y V2K - 1 

- < , < 1. 

z -^2K -I 

It must be that the curves never intersect. Thus, we can solve the Riemann problem 
< f/i, t/2 > in the rs— plane. 



43 

2.5. THE RIEMANN PROBLEM 




Figure 2.2. Possible Intersection of two 3— Shock Curves. 

Now, we use this result to find a solution to the Riemann problem with Ul = 
{Pl,vl, Sl) and Uji = {pr,vr, Sji). By the previous argument, find a middle state 
{Pm,vm) that solves the Riemann problem, < {pl,vl), {pl,vl) >, in the rs-plane. 
We only need to find the two values of S on either side of the contact discontinuity. 
This can be accomplished by determining the change in entropy, across the 1 and 
3— waves then adapting these changes to the left and right values of S. For example, 
the left middle state Um would have entropy value 5*^ if we had a 1— rarefaction 
wave, and would have entropy value Sm, where 5*^/ — Sl equals the corresponding 
increase in S across the shock wave. We can find the change in entropy by looking 
at the equation p/n"' = A{S) and solving for S. This is possible since A is strictly 
monotone increasing away from zero. Similar methods determine the value of S'j^ 
and the value of entropy in the right middle state. Since the entropy values of the 
middle states satisfy, Sl < Sm and Sr < S'j^j, we have Sm,S'j^^j > 0. The position 
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Figure 2.3. Solution to the Riemann Problem < Ul,Ur >. The 
states Um and t/^ differ only in S. 

of the entropy jump is determined by the particle path emanating from the initial 
discontinuity with speed vm- 

This construction determines the two unique states Um = {pm,vm, Sm) and 
U'm — {Pm, Vm, S'j^f) that solves the Riemann problem in the region p > 0, — 1 < < 
1 and ^ > 0. Figure EJl □ 

We parameterize the 1 — (resp. 3) shock/rarefaction curve by the change in r(resp. 
s) and define the strength of a shock or rarefaction wave as the difference in the values 
of either r for a 1— shock-rarefaction wave, or s for a 2— shock-rarefaction wave. We 
choose the orientation on our parametrization so that we have a positive parameter 
along the rarefaction curve and negative parameter along the shock curve. Therefore, 
the solution of the Riemann problem can be given as a sequence of three coordinates, 
(ei,e2,e3) where, ei denotes the change in the Riemann invariant r from Ul to Um, 
62 the change in 5* from Um to U'^j and 63 the change in the Riemann invariant s 
from U^j to Ur. In summary, for i = 1, 3 we have a shock wave of strength when 
ej < and a rarefaction wave of strength ej when > 0. 
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We adopt the following notation: 



a 
S 



Strength of 1 — Shock Wave 



Strength of 3 — Shock Wave 



Strength of 1 — Rarefaction Wave 



Strength of 3 — Rarefaction Wave 



Strength of Entropy S-Wave 



If (61,62,63) is the solution to the Riemann problem with states Ul, Ur, we would 
have: 



a 



61 < 
Otherwise 



/3 







£3 < 
Otherwise 



61 > 
Otherwise 



£3 




£3 > 
Otherwise 



We define 5 = E/j — E/, where E = ln(A(S')). The value of S may be recovered 
by recalling this definition and since E is a strictly increasing function of S by (lA3p . 
Also, we will denote 6^^ as the absolute change of E across a shock wave of strength u. 
More specifically, if two states were separated by a shock of strength u the absolute 
change in E across the shock would be S^j for either a 1 or 3— shock. Since we have 
shown that the change in E is independent on the base state and dependent only on 
the strength of the wave, 6^ is well defined. 
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6. Interaction Estimates 

In this section we prove estimates for elementary wave interactions with a method 
that follows the work by Nishida and Smoller, and Temple in [S] and [12] . This 
method is employed in order to simplify the estimates on the variation in the en- 
tropy. The alternative approach, useing the wave interaction potential Var {ln{p)} 
introduced by Liu and used in |10j , simplifies the estimates dealing with the first and 
third, nonlinear, characteristic classes, but complicates the estimates dealing with 
the entropy. 

Consider the following three states, Ul = {pl,vl,^l), Um = {pm,vm,^m), and 
Ur = {pr,vji, T,fj). We wish to estimate the difference in the solutions of the three 
Riemann problems < Ul,Um >, < Um,Ur >, and < Ul,Ur > with solutions 
denoted by a 1 subscript, 2 subscript and ' respectively. 

2.10. Proposition. Let Q be a simply connected compact set in rs— space. Then 
there exists a constant Cq, 1/2 < Cq < 1, such that for any interaction < Ul, Um > 
+ < Um, Ur >^< Ul, Ur > in Q at any value ofT,, one of the following holds: 

^.) A = < 0, 0<B< Co^, 

or 

5 = -e < 0, < A < CoC, 
a.) A<0, andB < 0. 
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Where A = a' — ai — 02 and B = (3' — (3i — (32 
and 3 shock waves in the solutions. 
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are change in the strengths of the 1 



Here we note that after an interaction, the shock wave strength in one family 
may increase, but this increase is uniformly bounded by a corresponding decrease in 
shock strength for the opposite family. 

Proof. These estimates are proven in Chapter [S] by a systematic look at all 
possible wave interactions. Because the interactions are independent of entropy 
level, we only consider interactions within the first and third characteristic classes. 
There are sixteen unique incoming wave configurations and between one and four 
possible outgoing wave configurations. The main idea is that after an interaction, 
there cannot he an overall increase in the strengths of the shock waves. This fact 
follows since as the solution progresses forward in time, cancelations and merging of 
shock and rarefaction waves of the same class lead to a decrease in shock strength. 
For example, when a shock wave is weakened by a rarefaction wave, a reflected 
shock wave is created in the opposite family. This interaction may increase the 
total strength of the shock waves in the opposite family, but the total gain in shock 
strength is uniformly bounded by the loss in the weakened or annihilated shock. 

We choose the constant Co to be the maximum slope of the largest shock wave 
curve that lies within the compact set i7 or 1/2 in order to bound the constant below. 
More specifically, let aJ be the strongest largest shock wave possible in VL. Then we 
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Figure 2.4. The creation of an entropy wave after elementary waves interact, 
take Co to be 



Finally, by Lemma 12.71 the slopes of the shock wave curves in a compact set in the 



For interactions in a compact set, the variation in S across a shock wave is 
uniformly bounded by a constant times the strength of the shock. But, the variation 
in E may increase after an interaction because of the likely creation of an entropy 
wave. Typically, across these waves the pressure is invariant and there is a jump 
in density; however, under the assumption fllUp . there must be no jump in energy 
density. Thus, we cannot use In (p/pi) or the change in the Riemann invariants r or 
s as a measure of wave strength. It should be noted that under certain interactions, 
such as an i— shock being weakened by an incoming i— rarefaction wave, an entropy 
wave is created with strength such that 5*^/2 — Smi is equal to the loss in entropy 



(46) 




rs— plane are strictly bounded away by 1. Therefore, we conclude Co < 1. 



□ 
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change across the shock, plus the change in the entropy across the new shock wave 
in the opposite family. We need a way to bound the variation in the entropy waves, 
and it turns out that this increase is bounded by a corresponding decrease in the 
shock strengths. 

2.11. Proposition. For every simply connected compact setQ inrs— space, there 
exists a constant M > such that after every interaction in fl, at any value S for 
the system ([6]) with ffTTl) in the ultra-relativistic limit, the following holds: 

\S'\ -\S,\- \S2\ + {S^, + 5^, - (5,0 + {5p, + - V) < -M{A + B). 

Proof. Choose Cq so that Proposition 12. lOl holds. Since is a compact set, let 

aJ = sup {||(ri, si) - (r2, S2)|| : (n, Si), (r2, S2) G Vt} . 

Then the strength of the largest shock wave in Vt is bounded by W. Furthermore, let 
M= (1 -Co)"^M, where 

(47) M = 2^E-^(^), 

duj 

which is twice the largest rate of change of S for all shocks contained in Vt. Also, 
since [S — Si](co') is positive and convex up, we have for strengths, oj' > Ui + c<j2. 

The proof will be split into two cases, one for each of the two cases from Propo- 
sition [2]T0l First let us assume that A <0 and B <0. i.e. 

a' — ai — a2 = —^a < and P' — Pi — P2 = — < 0. 
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We have, ai + 02 — = a' and hence, 6(^ai+a2-^a) = ^a'- It follows that 
Rearranging, 

(48) <5,, + 5^, - 5^, < ^M^a < -\mA, 
and similarly, 

(49) + 5^, - 5p, < hl^p < —MB. 

The right hand inequalities follow from the fact that M < M. Also, the change in 
entropy across the two Riemann problems before and the resulting one are equal: 



(50) Sa' + S' - 5fs' = 6a^ +61- 6fs^ + + 62-5^ 



/32- 



Rearranging (l50l) and using the previous estimates (HHjl and (l49l) . we find 

(51) {6' -61- 62) + (5^1 + Sp, - V) = « + S0.2 - So.') < -\mA 
and 

(52) (5' - 5i - 82) + (5,' - 5,, - 5aJ = (V - - 5^3,) > ^M5. 
Adding the inequality dH]) to flSTl) . 

(y - 5i - ^2) + (5a, + - 5«0 + {5p, + 5^, - V) < -^^^ - ^^^^ = -MA, 
and adding —1 times the inequality (H9ll to (l52l) . 

((5' - 5i - ^2) + (5a' - + 5aJ + (V - 5ft - 5ft) > ^MS + ^MS = MB. 
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Therefore, after multiplying the entire inequality by —1 we have 

- {S' -Si- 82) + + 5^, - <5„0 + (^/3i + -Sp,)< -MB. 

Since < —MA and < —MB by assumption, it follows that 

\6' -61- 62] + + 6^, - 5,0 + {6f3, + 5/32 - V) < -M{A + B). 

Furthermore, since \6' — 61 — 62] > — \6i\ — \62\, we deduce, 

-\6i\- \62\ + (da, + 6a, - 6a') + (6(3, + 6p, - V) < -MiA + B). 

This concludes the proof of the first case. 

Now, without loss of generality assume A = —^ < and < i? < Cq^. The 
mirror case when < A is similar. As before, we can obtain the estimates, 

6ai + 60,2 - 6a' < 

and 

(53) {6' -61- 62) + {6a, + 6a, - 6a') + (5/3, + 5/3, - 5^0 < Me 
From fl50|) we have, 

- (5' - 5i - 52) + (5,1 + 5„2 - 6a') + (5/3/ - 5/3, + 6^,) = 0, 

and since P' > Pi + P2, we have 5/3, + 6^2 — 5/3/ < 0, and so by adding this inequality 
twice, 

(54) - (5' - 5i - 52) + {6a, + 6a2 - 6a') + {6p, + 6^2 - 5^3') < 0. 
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Therefore, from (l53l) and fl5^ . 

which as before, reduces to, 

1^1 - l^ll - 1^21 + {Sa^ + Sa, - Sa') + {Sf}, + Sfi, - (5/3') < 

But, 

= M(l - Co)e = M(e - CoO < M{-A -B) = -M{A + 5), 

where we used the fact — Cq^ < —B following from the assumption that < -B 
and A = 
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CHAPTER 3 
The Glimm Difference Scheme 
1. Introduction 

In 1965, Glimm proved existence of solutions to general systems of strictly 
hyperbolic conservation laws with genuinely non-linear or linearly degenerate charac- 
teristic fields, [4]. To obtain existence, Glimm needed to restrict to initial data with 
sufficiently small total variation to avoid having to rule out the possibility that the 
complicated, global nonlinear structure of the conservation law might create finite 
time blow up of the solution or approximation scheme. His method takes a piece- 
wise constant approximate solution at one time step and uses numerous solutions to 
Riemann problems, defined at each point of discontinuity, to evolve the solution to a 
later time. After the approximate solution is brought forward in time, the solution is 
randomly sampled and a new piecewise constant approximate solution is obtained. 
A fascinating consequence is that one cannot choose any sequence of sample points 
to choose the states used for the new piecewise constant function at each time step, 
but rather must sample outside a set of measure zero in the space of all possible 
choices. 

One way this scheme may break down for general systems of hyperbolic conserva- 
tion laws is that Riemann problems may not have solutions if the initial left and right 
states are sufficiently far apart. A canonical example of this phenomenon occurs in 
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the p— system which models a classical isentropic gas in Lagrangian coordinates. For 
this system, if the difference in velocity of the two initial states is sufficiently large, 
all the gas will be pulled from the region in between the two states forming a vacuum, 
[9] . This possible complication and issues with large scale non-linearities, led Glimm 
to prove existence for initial data with small variation. He showed that in this case, 
the total possible increase in variation in the approximate solution is bounded by 
a corresponding decrease in a quadratic functional. Thus, having a bound on the 
total variation in the solution showed that the Riemann problems used in evolving 
the approximate solution can be defined for all time. 

In our case we prove a large data existence theorem; there is no restriction on the 
"smallness" of the initial conditions. In our existence proof, we will not need to use a 
quadratic functional to bound the total variation, because the geometric structure of 
the shock and rarefaction curves in rs— space do not allow the approximate solution 
to behave badly in the large. In this section, we will introduce the Glimm scheme 
and use it to construct solutions to ([6]). 

2. Glimm Difference Scheme 

We say U{x,t) is a weak solution of with initial data Uo{x), if for all ip G 
[M+, M] the following holds: 
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■X- 



(fc + 2)Ar 



(A; + l)Zk 



Figure 3.1. Construction of Piecewise Constant Function 



We begin by partitioning space into intervals of length Aa; and time into intervals 
of length At. In order to keep neighboring Riemann problems from colliding, we 
impose the following CFL condition: 



Note for 1 < 7 < 2 this condition is satisfied since the characteristic speeds IHM are 
bounded above and below by 1 and —1. 

We inductively define our approximate solution. To begin suppose that we have 
an approximate solution at time t = nAt, U{x,nAt), which is constant on the 
intervals, {kAx, {k + 2)Aa;), where k + n is odd. At each point x = kAx a Riemann 
problem is defined. Solve each Riemann problem for time t = At. This evolves our 
approximate solution from t = nAt to t = (n + l)At. To finish, we must construct 
a new piecewise constant function at time t = {n + l)At. Choose a G [—1, 1] and 
define, U{x, {n + l)At) = U{{k + 1 + a)Ax, {n + l)At-) for x G {kAx, {k + 2)Ax) 
and k + n + 1 odd. The term At— denotes the lower limit. 




I = 



1,2,3. 




(A;-2)Z^ (fc-l)Zk AAr {k+2)/^ 

Figure 3.2. I-Curve J. 



To begin this process ai t — 0, obtain a piecewise constant function from the 
initial data Uo{x) by again choosing a e [—1, 1] and defining, U {x, 0) = Uo{{k+a)Ax) 
for k odd. 

Consider, 9 e Hi^o [~^' ■'■^ other words, 6* = (^O) ^i, ■ ■ ■ ,&n, ■ ■ •) with 9i e 
[— 1, 1]. Then, for initial data, we say Ue^Axi^^t) is the approximate solution given 
by a mesh size of Ax with sampling points at the n}^ time step given by 9n- 

In order to estimate the change in the variation of our approximate solutions, we 
will define piecewise linear, space-like curves, called I-curves, which connect sample 
points at different time levels. If an I-curve J passes through the sampling point 
{{k + 9n)Ax,nAt), then on the right J is only allowed to connect to ((A; + 1 + 
9n±i)Ax, (n ± l)At) and on the left to {{k - 1 + 9n±i)Ax, {n ± l)At). 

We consider two functionals defined on 7— curves and will analyze how the func- 
tional change as we change from one /—curve to another. This will allow us to 
estimate the change in variation of the approximate solution as it is evolved using 
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the Glimm scheme. We define for an J— curve J: 

(55) F(J) = 5^a, + 5^A + V^ 

J J 

and 



(56) L( J) = 5^ (a, - Mo5„J + (A - Mq^^) - Mq |5| + \/, 

J J j 

where the sums are taken over all waves, or fractions of them in the case of rarefaction 
waves, that cross J. The constant Mq will be chosen later and V = Var {Uo{-)} is 
the variation of the initial data. 

The main problem in our analysis is to show that the variation in the entropy 
waves stays bounded for all time. To do this we need to bound the possible change 
in E across shock waves. This is accomplished by first showing that the variation in 
r and s stays finite for all time. This implies that all the interactions, as projected 
onto the rs— plane, occur in a compact set. Thus, there is a largest possible shock 
strength in this compact set, and using the fact that the derivative of the entropy 
change as a function of wave strength is monotone increasing, there is a constant 
such that the entropy change is bounded by a constant times the wave strength. 
Moreover, we can then use Proposition 12. Ill to estimate the increase in the variation 
in entropy in our approximate solutions. 
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3. Estimates on Approximate Solutions 

For initial data Uo{x) and corresponding approximate solution Ug^Ax{x,t), define 
Uo{x) and Uo,Ax{x,t) as the initial data and approximate solutions viewed as func- 
tions of r and s only. The first estimate will show that the variation in the Riemann 
invariants across an I-curve J is bounded above by the functional F(-) on J. 

3.1. Proposition. Let Uo{-) be of finite variation. If the approximate solution 
Ue,Axix,t) is defined on an I-curve J, then, 

(57) Varrs{J) < ^F{J). 

Proof. Let Var~{J) denote the variation across J given by a decrease in r. The 
only waves that contribute to the decrease in r are 1 and 3— shocks. Furthermore, 
we have 

(58) l^ar;(J)<^a, + 

J J 

where the sum is over all waves of the particular type crossing J. We can similarly 
define Var^{J) as the variation given by increases of r across elementary waves. The 
only increase is given by 1— rarefaction waves, 

(59) yar+(J) = ^//i. 

J 

Following this line of reasoning for s, we also have 

(60) Var:{J)<Y,Oi^ + Y,l^i 

J J 



59 

§3.3. ESTIMATES ON APPROXIMATE SOLUTIONS 
and 

(61) Vart{J) = Y,Vi- 

J 

The initial data Uq may be written as a function of the Riemann invariants r and 
s, Uq{x) = (ro(x), so(x)). Since f/o(') is of finite variation, the following limits must 
exist: 

lim ro(a;) = r^, lim so{x) = s^. 

X— >itoo a:— >itoo 

Indeed, let {x„}^q be an increasing sequence of real numbers such that ^ oo as 
n — > oo. Then, 

oo 

\roixn) - ro(x„_i)| < Var{Uo{-)} < oo. 

n=l 

Hence the sequence ro(a;„) is Cauchy, which converges to a finite limit The other 
cases are entirely similar. 

For any I-curve J, the end states at ±oo are given by (r^,s^). From this we 
obtain 

\Var+{J) - Var;{J)\ = |r+ - r"] < V, 

and hence, 

Var+{J) < Var;{J) + V. 

Using (1581) and (|59]). 

J J J 

Similarly from ( l60i) and ( l6Ti) . 

J J J 
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Combining these together we have, 

J 

Thus, 

VaVrsiJ) 



□ 

We will now show that the functional F(-) on the I-curves is non-increasing. We 
define a partial ordering on the I-curves by saying that J -< J' if the curve J' never 
lies below the curve J. Furthermore, we say that J' is an immediate successor to J if 
J ^ J' and J and J' share all the same sample points except for one. It is clear that 
for any pair of I-curves such that J ~< J', there is a sequence of immediate successors 
that begins at J and ends at J'. The next proposition shows that if our approximate 
solution is defined on an I-curve, it can be defined for all following I-curves. 

3.2. Proposition. Let J and J' be I— curves, J -< J' , and suppose that J is 
in the domain of definition of Ug^Ax- If F{J) < OO; then J' is in the domain of 



IJ'i 



J \ J J / 



\ J J J J J 



< 4:F{J). 




(/c-l)Zk hAc (A; + l)Zk (k + 2)Ai: 



Figure 3.3. Immediate Successor I-Curves J' and J. 

definition ofU/^xfi, and F{J') < F{J). Moreover, if Var^s {Uo{-)} < oo then t/^^Ax 
can be defined for t > 0. 



Proof. We proceed by induction. Suppose first that J' is an immediate succes- 
sor to J. Then the difference F{J') — F{J), is given by the change in shock wave 
strengths across the diamond enclosed by J' and J. See Figure 13.31 This is a conse- 
quence of the fact that the waves the head into the diamond from the left and right 
solve the same Riemann problem as the outgoing waves in the new single Riemann 
problem. If we denote Jq and Jo as the diamond portion of J' and J, we have, 

F{f)-F{J) = + + + J]/3, + V^ j , 

J' J' \ J J J 

Jq Jq Jo Jo 

= (a -ai- tta) + iP' - A - P2) , 
= A + B <0. 
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The last line follows from Proposition 12.101 Thus, F{J') < F{J) for immediate 
successors. For any a general J and J' such that J -< J', we produce a sequence 
of immediate successors that take J to J'. At each step the functional F is non- 
increasing, thus F{J') < F{J) continues to hold. 

By Proposition 13.11 Varrs{J') < ^F{J') < 4:F{J), so, J' is in the domain of 
definition of Ue,/\x- Moreover, if Vavrs {f^o(')} < ^5 then Varrs{0) < 00 for the 
unique /—curve that lies along the line t = 0. In order to show that U/\x,e can 
be defined for t > 0, we must show that VaVrs {f^6»,Ax(-, i)} < 00 for all time. But, 
this condition is equivalent to showing the variation across any /—curve J is always 
finite. Since for any /—curve J, 

VaVrsiJ) < ^F{J) < 4F(0) < SVar^s {Uo{-)} , 

the result follows. □ 

Again, Proposition 13.21 shows that the variation of our approximate solution in 
the variables r and s is finite. Thus, there exists a compact set in the rs— plane that 
contains all the interactions in our approximate solution. 

3.3. Corollary. Suppose that Vavrs {Uq{-)} < 00. Then there exists a sim- 
ply connected compact set Q in the rs— plane such that all possible interactions are 
contained in Q. 

Proof. From Proposition 13.11 and Proposition 13.21 we know that for any I-curve 

J, 

VaVrsiJ) < ^F{J) < 4F(0) < 8V^ar,, {Uo{-)} = N < 00. 
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Thus, the distance between any to states occurring anywhere in our approximate 
solution is bounded by A^. Consider the left limit state of f/o(')> ('^ 5 )■ Therefore, 
all states must be contained within, i?2Ar(?"~, the ball of radius 2N centered 
around (r^, s^). □ 

Now, we show that the variation of our approximate solution, including the 
variation in E, is bounded above by the functional L{-). 

3.4. Proposition. Suppose Var {Uo{-)} < 00 and J is an I— curve that is in 
the domain of definition of Ug^Ax- Then there exists constants Mq > and K > 0, 
independent of Ax and 6, such that, 

(62) Var{J) <K ■ L{J). 

Proof. The variation across the /—curve J is bounded by 

Var{J) < \^ar(Shock Waves) + Var (Rarefaction Waves) 
+V^ar(S- Waves) + VariTj across Shocks). 

Since VaVrs {Uq{-)^ < Var {f/o(-)} < 00, we have from Corollary 13.31 that all the 
interactions projected into the rs— plane occur in a compact set Vl. Therefore there 
exists a constant M > such that for a shock wave of strength cj, 5^) < Mu. Let 
M = (1 — Cq)~^M as in Proposition 12.111 Since, M < M we have for a shock wave 
of strength uj, 5^ < Moj. 

From the proof of Proposition 13.11 we can bound the variation from the shock 
waves and rarefaction waves by the shock waves crossing J and the initial variation 



§3.3. ESTIMATES ON APPROXIMATE SOLUTIONS 
V. Thus, 

Var{J) < 2 [^a. + ^A] + 

\ J J J J J 



\j J J J \ J J / 

< (4 + M) ( + + + 

\ J J J J 



Let Mo < 1/2M. Then, 



Thus, for a shock wave of strength a;, 

u;<2{uj- Mq5^). 

Using this, we find, 

Var(J) < 2(4 + M) ( (a, - Mo5„,) + (A - Mo^^^J + ] + |5| 
\ J J / J 

Finally, we put the sum of the strengths of the entropy waves inside, 

Var{J) < 2(4 + M) ( 5^ (a, - M^Sp,) + (A - Mq^^J + Mq |<5| + 
\ J J J 

We can do this because. 

Mo • 2(4 + M) > 2MMo > 1. 
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Therefore, 



with K = 2(4 + M). 



Var{J) <K-L{J), 



□ 



3.5. Proposition. Suppose that Var {Uo{-)} < oo and J, J' are I-curves such 
that J -< J' and L{J) < oo. Then J' is in the domain of definition of U0^/^x{x,t) , 
L{J') < L{J) and UAx,e{x,t) is defined for t > 0. 

Proof. Since Varrs {Uo{-)} < Var{Uo{-)} < oo there exists a compact set O 
that contains all possible interactions. Define M as in Proposition 12.111 and take 
Mo < 1/2M. As with Proposition 13.21 we prove the result by induction on the / 
curves. First let J' be an immediate successor to J. Let Jq and Jo be the parts of J' 
and J that bound the diamond formed by J and J'. Using this and the definition 
of L(J), 



L{J')-L{J) < 



4 



4 



J2 («. - Mo6^J + 5^ (A - Mo6fsJ + Mo 5^ \6\ 

. Jo Jo Jo 

+Mo {6f3, + - V) + Mo - \5i\ - \52\) . 



Now we refer to Proposition 12.101 and 12.111 We see that the first two terms are equal 
to {A + B) and the others are bounded above by —M{A + B). Putting this together, 

L(J') - L{J) <{A + B)- MMo{A + S) < ^(A + 5) < 0. 
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For immediate successors, we have L{J') < L{J). Moreover, by Proposition 13.41 we 
have that the variation along J' is bounded by L{J') and hence L{J). Thus, J' is in 
the domain of definition of U/^x,9- 

For general J and J' such that J -< J', the same conclusion holds by constructing 
a sequence of immediate successors to move from J to J'. Along each step, the results 
above continue to hold. 

Finally, if Var {Uo{-)} < oo, we have L(0) < oo and for any /—curve J, L{J) < 
L(0). Which we can conclude that 

Var{J) < 2(4 + M)L{J) < 2(4 + M)L(O) < oo, 

so our approximate solution can be defined for t > 0. □ 
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CHAPTER 4 
Existence Theorem for Two Gasses 

In this chapter, we use Ghmm's Theorem [4] to prove existence of solutions 
to ([6]) in the ultra-relativistic hmit with an equation of state of the form (JTTl) . It 
should be noted that for 6 fixed and x„ = 1/2", the set of approximate solutions 
{f^e,Aa;„(3;, t)}^^ has uniformly bounded variation by Proposition 13.41 Furthermore, 
since the variation is bounded and each approximate solution has the same limits at 
infinity, the sup norm is also uniformly bounded. The approximate solutions are 
Lipschitz in time too since 

\\Ue,/ix{,-,t) — Ue^/s^xi,-, s)\\-i,i < C {Sum of all wave strengths} ■ 

{Maximum Speed of Wavefronts} < C'\t — s\. 

At this point Helly's Theorem ^ provides a convergent subsequence, t/e,Aa:„- (a;, t), 
that converges to a function U{x, t) with finite variation for each fixed time. However 
at this time, there is no justification that this limit function is actually a weak so- 
lution. Glimm's Theorem guarantees that there exists a subsequence that converges 
to a weak solution. 
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1. Existence of Weak Solutions 

3. Theorem (Glimm, 1965). Assume that the approximate solution Ue^^xi sat- 
isfies, 

(63) Var{Ue,AxX-,t)} <N <oo 

for Xi = 1/2^ , 6 E A = YYi^Q, and all t > 0. Then there exits a subsequence of mesh 
lengths Axj^^ such that Ug^/\xi^, ^ U in where U{x, t) satisfies, 

Var{U{-,t)} < N. 

Furthermore, there exits a set of measure zero A C A such that if 9 E A — A then 
U{x,t) is a weak solution to IQ. 

We now prove Theorem [T] by showing that our approximate solutions meet the 
assumptions of Glimm's Theorem. 

Proof. Assume the initial data satisfies, dHD, (fT5|) . and (fT6i) We will show that 
for all Axi and sample points 6, 

(64) Var{UA.A-,t)} <N <oo, 

where Ug^Ax{p{x,t),v{x,t), S{x,t)) = {Ui,U2,U^)0^ax- First we show that the varia- 
tion in p, f , and S is bounded for all time in the approximate solutions. 

From Proposition 13.11 and Proposition 13.21 we have that the variation of our 
approximate solution in r and s is uniformly bounded for all time. More specifically, 
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Varrs{Ue,Ax{-:t)) < 4F(0), 



< 4 



^Q;i + ^A + yar,,(C/o) 



< S-Varrs{Uo{-)). 



Prom this estimate, we show that the variation of 



In (p) and In 



1 + v 
1-v 



are also bounded for all time. Indeed by the definition of r and s, 



we have 



Var < In 



1-V{;t) 



]-supy^ \{r{xi+i,t) + s{xi+i,t)) - {r{xi,t) + s{xi,t))\, 
1 ^ 

1 

1 ^ 

+ ^sup^ - s{xi,t)\, 



N 



1=1 



< 8-Var{Uo{-)} 



Similarly, using 



ln(p) = ^-t^(s-r). 
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we find, 



Var{Hp{;t))} < 16 ( ) Var{Uo{-)} 



Now, we show the variation in E is bounded for all time in approximate solutions. 
This is clear from Proposition 13. 41 and Proposition 13. 51 because there exists a constant 
M so that 

Var {Ee,Ax(-, t)} < 2(4 + M)L(O). 

We can now show that the variation in p, v and S is bounded for all time. Since 
Var {ln(p(-, t))} < oo for alH > there exists a constant 6 > such that p{x, t) < b. 
Let c = max {!,&}, then 

N 



Var{p{-,t)} = snpy^ \p{xi+i,t) - p{xi,t)\, 

N ^ 



i=l 



N 



< c- sup I ln(p(xi+i,t)) - ln{p{xi,t))\, 

N ^ 

1=1 

< c-Var{ln{p{-,t))} . 



For V we have, 



TV 



Var{v{-,t)} = snpy^ \v{xi+i,t) - v{xi,t)\ , 

N ^ 

1=1 

i=l ^ 



1, 



< -Var < In , , , 

2 1 \l-v{-,t) 



v{Xi+i,t) 

l + v{-,t) 



-In 



1 + v{xi,t) 

1 - v{Xi,t) 



The factor 1/2 comes from the fact that the slope of the chord connecting the points, 



v{xi+i,t),\n 



1 + v{Xi+i,t) 
1 - v{Xi+i,t) 



and ( v{xi, t), In 



1 + v{xi,t) 

1 - v{Xi,t) 
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is bounded below by 2. 



For S we need to find a constant C such that 



\S{x,t)-S{y,t)\<C\^{x,t)-^iy,t)\. 



Since S is of finite variation for all time, there exists a largest and smallest value of 
S, say Smax and Smin with < Smin < Smax- Define C by 



Finally, from Proposition 12.51 the determinant of the Jacobian is bounded away 
from zero for all approximate solutions. Therefore, the variation in conserved vari- 
ables, (?7i, f/2, U3), are bounded for all t > 0, 9 and Axi. 

Therefore, Theorem [3] provides existence of a set measure zero A C A such that 
if we choose 6 E A — A there exists a subsequence of mesh refinements, Axj^. 
such that ?76»,A^. converges pointwise almost everywhere in L]^^ to a weak solution, 
U{x,t) of IQ. Moreover, this solution satisfies 




It follows that. 



Var{S{-,t)} < C -Var t)} . 



Var {ln(p(-,t))} < N, 




and 



Var{S{-,t)} < N, 



§4.1. EXISTENCE OF WEAK SOLUTIONS 
for some > 0, alH > and is Lipschitz in 



73 



CHAPTER 5 
Interaction Estimates 

We give a systematic approach to the wave interaction estimates needed to prove 
Proposition 12.101 From the special geometry of the shock-rarefaction curves in the 
space of Riemann invariants we can analyze the interactions as done for the classical 
Euler equations in |8] and [12j. There are sixteen possible incoming wave profiles, 
and among these one to four different outgoing wave profiles, each of which will be 
covered on a case by case basis. We assume that all the interactions occur in a simply 
connected compact set f2 C M^. Recall that for a compact set fl we have defined the 
constant Cq, (H6|) . as the max of 1/2 or the largest slope possible of a shock curve 
contained in Q. Since the shock wave slopes are strictly bounded above by 1, we 
have, < 1/2 < Co < 1. 

During these estimates we repeatedly use the fact that the shock curves in the 
space of Riemann invariants are translationally invariant, convex and whose deriva- 
tives are bounded above by a constant. We reference Lemma [2.71 for these results. 
Since our definition of wave strength is determined by the change in r for 1— waves 
and s for 3— waves, we use the following two facts: 




Figure 5.1. Shock Curve Slopes are bounded by Cq. 

(1) The change in s along a 1— shock is uniformly bounded by the change in r 
and vice versa for 3— shocks. Indeed, since we have 



ds ^ V2K - 1 



dr -J2K - 1 



for 1— shocks and 



dr V2K- 1 

< 



ds -V2K-1 
for 3— shocks, we have for our constant Cq, 



— < Co and — < Cq. 

Zl Z3 



See Figure [STTl 

(2) Suppose two shock curves of the same family that begin at two distinct 
states Ui and U2 and meet at a common third state U3. Then the ratio of 
the distances along the r and s axes from Ui and U2 are bounded by Cq. 




Figure 5.3. + (02, /?2) ^ A<0, B>0 

Again, we have, 

— < Co and — < Cq. 

See Figure 15.21 
We now begin our interaction analysis. 

(1) + (a2,/?2) 

• If A < and i? < we are done. 
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Figure 5.4. + (02, /?2) ^ {a',P'), A>0, B <0 



+ {a2,f32) ^ Figure 
Suppose that A = < and B >0. We have A = yi — ^2 and 
B = 1/2 — zi. From A < 0, ?/2 > -^i, and hence, 



yi < Zi < 1/2 < Z2. 



Therefore, 



A + B = {yi-z2) + {y2-z,), 



{yi + y2) - (^1 + 2:2), 



< (Co-l)(Zl+^2), 



< (Co-l)(^2-Z/l), 



< C^i + A. 



Hence, B < CqC,. Note, the inequahties hold since (Co — 1) < 0. 




Figure 5.6. + (02,^2) {a',v') 

A) + (a2,/92) ^ {a',P'y. Figure [El 

Suppose that A > and -B = — ^ < 0. We have a' = ai + 
Hi + a2 + 1/2 and hence, A = a' — a 1 — 0:2 = yi + y2- Furthermore, 
l3' + Z2 + zi = l3i + (32, which gives, B = (3' - jSi - (32 = -zi - Z2 = 
Thus, A < Co^. 

+ (02, /?2) ^ Figure [53] 

A = -z, 5 = y < CqZ. 
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V2 



Figure 5.7. + (02, ^72) ^ 




Figure 5.8. + (a2,?72) 

• + (a2,/52) ^ (a', 77'): Figure [5J] 

B = —Pi — P2) = —z, and yl = a' — cti — 02 = ?/ < Cqz. 
(2) (ai,/5i) + (a2,r/2) 

• (ai,/3i) + (a2,?72) ("', Figure [521 

B = —z, A = a' — ai — a2 = y < Cqz. 

• (ai, A) + (a2,?/2) ^ («',/?'): Figure ESI 

= P'—Pi—P2 = —-2- Since, ai+a' = a2+y+ai =^ a'—a2 = y 
then, A = a' — ai — a2 < y < Cqz. 





(3) (ai,/5i) + (/i2,/52) 

• + (/i2,^2) ^ Figure ED 

A = -ai = -z, B = p' - P^- P2 = y < CqZ. 

• + (/i2,/32) Figure EIO] 

A = a' - ai = -z, B = - I3i - 132 = y < Cqz. 

(4) (ai,/?i) + (/X2,r/2) 
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a/ 


m 






y 


M2 



z 



Figure 5.11. + (/i2,^?2) 




Figure 5.12. + (/i2,r?2) ^ (/"',/?') 

• (tti, A) + (/i2,r/2) 

A < 0, 5 < 0. 

• (ai,/3i) + (/i2,?72) ^ Figure |MU 

A, 5 < or, 5 = = -2: and A = a' - ai<y < Cqz. 

• (ai, A) + (^2,r?2) Figure [5l2 

A,B < or; A = -z and y + /?i = /?' + 772 ^ - /?i = ?/ - 772 
^ 5 = /3' - A < y < Co^. 
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Figure 5.13. (ai, A) + (/i2,^?2) {a',l3'), < A and < fi. 




Figure 5.14. {ai,r]i) + {a2,P2) («',/?') 



• (ai, A) + {l22,V2) («',/?'): Figure [513] 

B = P' — Pi = —z and ai + ?/ = /i2 + a' 
A = a' — ai < y < Cqz. 

• + (/i2,r/2) Figure [5l3] 

A = a' — «! = —2; and y + Pi = P' + Tj2 ^ P' — Pi = y — rj2 
B = P'-Pi<y< Cqz. 
(5) {ai,r]i) + (02, /32) 




Figure 5.15. (ai,?7i) + (02,^2) 




ai 



Figure 5.16. (ai,?7i) + (a2,?72) («', V) 

• {ai,T]i) + (02, /?2) ^ Figure EH 

A = a' — ai — ^2 = y and i? = — /32 = — z ^ A < Cqz. 

• {ai,r]i) + {a2,(32) — > {a',v')- Figure [5J5] 

A = a' — ai — 0:2 = y and B = —(32 = —z =^ A < CqZ. 

(6) (ai,?7i) + (a2,?72) ^ ia',r]'): Figure |5J6] 

A = a' - ai - a2 = < and 5 = < 

(7) (ai,r?i) + (/X2,/32) 



/32 



m 



/i 



ai 



Figure 5.17. (ai,r/i) + {^12,^2) (yu',^') 




7/ 




Figure 5.18. {ai^rji) + (/i2,/?2) ^ ^,5 < and < A. 



(ai,?7i) + (/i2,/32) (yu',??'): Figure |5J7| 
A = -ai < and = -/^s < 0. 

(ai,?7i) + (/i2,/32) (a',??'): Figure |5J8| 

Either A < and 5 < or A < y and B 

A < -CoB. 



-P2 = z, therefore, 




Figure 5.19. {ai,rii) + {f^2,p2) {ft',P'), A,B <0, andO < B 



Figure 5.20. (ai,?7i) + {^12,^2) Case 1. 

• (ai,77i) + (/i2,/?2) ^ Figure Em 

Either A < and B < or B = (3' -(32 <yandA = -ai 
therefore, B < Cqz. 

• (ai,r/i) + (/i2,/92) ^ Figure EID and ESQ] 

Case 1.) A<0 and B <0. 

Case 2.) A = a' - ai = -z and B = (3' - (3i < y < Cqz. 
Case 3.) B = P' — P2 = —z and A = a' — ai = y < Cqz. 





Figure 5.21. {ai,r]i) + {n2,P2) {a',P'), Case 2 and 3. 




z 



Figure 5.22. (ai,?7i) + (/i2,^?2) ^ (/i',/?') 
(8) {ai,r]i) + {iJ,2,V2) 

A = -ai < and 5 = < 0. 

^ = -/i2 < and 5 = < 0. 

• + (/i2,r/2) ^ Figure [521 

A = — «! = —2; and B = (3' < y < Cqz. 

• (ai,?7i) + (/i2,r?2) ^ (a',/?'): Figure [523] 

A = a' — ai = —z and B = [3' <y < Cqz. 




z 



Figure 5.23. + (/i2,^?2) 

[J Ml 




Figure 5.24. + (02, /32) ^ (/i',/?') 

(9) (/ii,/5i) + (a2,/52) 

• + (a2,/32) Figure [521 
y4 = — 0^2 = —z and B = y < Cqz. 

• + (a2,/52) ^ Figure [E25] 

y4 = a' — 0^2 = —2; and B = (3' — [3i — [32 = y < CqZ. 
(10) (/ii, A) + (/i2,/32) ^ (/i',/?'): Figure [52a 

A = < and 5 = - /?i - /?2 = < 0. 



Figure 5.25. + (02,^2) {a',P') 



III 




fh 








P2 





Figure 5.26. (fiuPi) + (112,^2) ^ 

(11) + (02,772) 

A = -as < and = < 0. 
• (/ii,/?i) + («2,r/2) ^ (/i',/5'): Figure ESZl 

A, 5 < or, A = -a2 = -z and B = f]' - f3i < y < Cqz. 




Figure 5.28. (/ii, A) + (02,^72) {a',r]') 

(/ii, A) + {a2,V2) Figure [E28] 

A,B<{] or, we have y + 02 = A^i + Therefore, with i? 
— /?! = —2;, A = a' — a2 = y — fii<y< Cqz. 

+ (^2,772) ^ Figure [E29] 

Case 1.) A < and 5 < 0. 

Case 2.) A = a' — ^2 = —-2- Then we know [3' + rj2 = l3' + y a. 
therefore, B = p' - f3i < y < Cqz. 




Case 3.) B = (3' — (3i = —z. Then we know a' + fii = 0I2 + V 
and therefore, A = a' — 02 < y < C^z. 
(12) (/ii,/5i) + (/i2,r/2) 

• + (/i2,??2) ^ (yu',??') 

A = < and 5 = < 0. 

• + (/i2,r72) ^ (/i',/?'): Figure EJO 
A = and 5 = - /52 < 0. 
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Figure 5.31. + (/i2,?72) («',??') 




M2 2/ 



Figure 5.32. + (/i2,^?2) ^ («',/?') 

• + (Ai2,?72) ^ (a',??'): Figure ESU 

B = -Pi = -z and A = a' < y < CqZ. 

• (/ii, A) + (/i2,r72) ^ Figure Eja 

We have B = P' — Pi = —z and y = a' + /ii + /i2. Thus, 
A = a' <y < CqZ. 
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"2 




V2 



Figure 5.33. {fii,r]i) + (02,^2) (yu',??') 




Figure 5.34. {fii,r]i) + (02,^2) {a',r]') 

(13) (/ii,?7i) + ia2,f32) 

• (/ii,?7i) + {a2,P2) Figure [5J3] 

We have A = -02 < and 5 = -p2 < 0. 

• (Aii,?7i) + (a2,/32) {a',r]'): Figure [SJl 

i? = — = —z < and y4 = a' — 02 = ?/ — /ii <y< CqZ. 

• {fii,r]i) + (a2,/32) ^ (/i',/?'): Figure [OS 

A = -as < 0, 5 = /?' - /?2 < -^/2 < 0. 



"2 





(^2 















Figure 5.35. + (a2,/32) ^ 




^1 



Figure 5.36. {\i\,ri\) + (02, /92) (a',/?') Case 1 and 2. 



• (/ii,?7i) + (a2,/52) {ci\(3'): Figure [536] 

Case 1.) B = P' — P2 = —z. Since 0:2 + 2/ = /^i + a', we have 
A = a' — a2 = y — fJ'i- Therefore, A = a' — a2 < y < CqZ. 



Case 2.) A = a' — 02 = — 2- Since P2 + y = P' + ^1, we have 
B = /3'-l32<y< Coz. 



M2 



m 




! 








III 



















(32 



Figure 5.37. (/xi,r/i) + (/i2,/?2) ^ (z^',^') 



M2 



m 






Ml 


m' 




1 


» 1 

-^z / / 



Figure 5.38. (/ii,?7i) + (1^2,^2) (/u',/?') 

(14) (/ii,r/i) + (/i2,/52) 

• (Aii,?7i) + (/i2,/52) il^'^v')- Figure [533 

We have A = < and 5 = -/^a < 0. 

• (/ii,r/i) + (/i2,/32) {fi',P'y. Figure [ESa 

A = < and 5 = /3' - /32 = -r/i < 0. 

• (/ii,r7i) + (/i2,/?2) ^ {(y',v')- Figure [539] 

-B = — = —-2 and A = a' < y < Cqz. 





m 












y 













z 



Figure 5.39. (Aii,?7i) + (/i2,/32) ^ (a',??') 




Figure 5.40. (/xi,r/i) + (/x2,/52) 

• + (/i2,/?2) ^ (a',/5'): Figure [52Q] 

5 = - /?2 = and A = a' < y < C^z. 

(15) (/il,?7l) + (/U2, 772) 

• + (/i2,^?2) ^ (/i',??'): Figure [531] 
We have A = < and 5 = < 0. 
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V2 








m 



















Figure 5.41. (/xi,r7i) + (/X2,^2) ^ 





<^2^ 


/i' J 




Ml 



Figure 5.42. (/ii,r/i) + {a2,m) {l^'.v') 




VI 



Figure 5.43. {fii,r]i) + {a2,m) ^ 

(16) (/ii,?7i) + (a2,'72) 

• (Aii,?7i) + (a2,?72) {f^',v')- Figure [532] 

We have A = -ai < and 5 = < 0. 

• (/ii,r7i) + (a2,^?2) {a\v')- Figure [533] 

We have A = a' - a2 = -/ii < and 5 = < 0. 



m 



Figure 5.44. (/ii,r/i) + (02,^2) (yu',/?') 




Figure 5.45. + (02,^2) {a',l3') 

(/xi,r/i) + (02, V2) {^J-\l3'y. Figure [EH 

We have A = —a2 = —z and B < y < Cqz. 
(/ii,?7i) + (a2,?72) Figure |M5| 

We have A = a' — a2 = —z and B = (3' < y < Cqz. 
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